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A hydrophobic constriction site can act as an efficient barrier to ion and water 
permeation if its diameter is less than the diameter of an ion's first hydration 
shell. This hydrophobic gating mechanism is thought to operate in a number 
of ion channels, e.g. the nicotinic receptor, bacterial mechanosensitive channels 
(MscL and MscS) and perhaps in some potassium channels (e.g. KcsA, MthK, 
and KvAP). Simplified pore models allow one to investigate the primary charac- 
teristics of a conduction pathway, namely its geometry (shape, pore length, and 
radius), the chemical character of the pore wall surface, and its local flexibility 
and surface roughness. Our extended (ca. 0.1 fis) molecular dynamic simula- 
tions show that a short hydrophobic pore is closed to water for radii smaller 
than 0.45 nm. By increasing the polarity of the pore wall (and thus reducing 
its hydrophobicity) the transition radius can be decreased until for hydrophilic 
pores liquid water is stable down to a radius comparable to a water molecule's ra- 
dius. Ions behave similarly but the transition from conducting to non-conducting 
pores is even steeper and occurs at a radius of 0.65 nm for hydrophobic pores. 
The presence of water vapour in a constriction zone indicates a barrier for ion 
permeation. A thermodynamic model can explain the behaviour of water in 
nanopores in terms of the surface tensions, which leads to a simple measure 
of "hydrophobicity" in this context. Furthermore, increased local flexibility de- 
creases the permeability of polar species. An increase in temperature has the 
same effect, and we hypothesise that both effects can be explained by a decrease 
in the effective solvent-surface attraction which in turn leads to an increase in 
the solvent-wall surface free energy. 
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ohver@biop.ox.ac.uk, WWW: |http : //san som . biop . o x . ac . uk/j 



1 



1 Introduction 

Although only about 3 nm thick, the mem- 
branes of living cells present an efficient 
barrier to polar substances such as wa- 
ter or ions. Transport of water and ions 
into and out of the cell is facilitated by 
specialised proteins. For molecules that 
move down their electro-osmotic gradients 
the proteins are ion channels 0] and wa- 
ter pores (aquaporins) 0. They share 
common functional characteristics such as 
high specificity (for instance, the potas- 
sium channel KcsA prefers potassium ions 
over sodium ions, aquaporins allow the 
flow of water but not of ions or protons) 
and transport rates comparable to diffu- 
sion in bulk solution. Ion channels can also 
be switched ("gated") between an open 
(ion conducting) and closed (ion blocking) 
state by external signals such as changes in 
transmembrane voltage, binding of a lig- 
and, mechanical stress etc. 

In recent years a number of near atomic 
resolution structures of ion channels 0- 
[13] and aquaporins 11- 13] have been pub- 
lished. Based on these structures an atom- 
istic understanding of the transport and 
gating properties is emerging. If we view 
transport proteins as molecular machines 
which are designed by evolution to perform 
selective and efficient transport that can 
be controlled, i.e. gated, then we can ask 
what the "building blocks" of these ma- 
chines are. This involves not only iden- 
tification and structural characterisation 
of the protein domains involved in gat- 
ing, but also, at a more abstract level, un- 
derstanding the underlying physical prin- 
ciples. 

Here we address the questions of what 
are physical properties that are important 
for the permeation of ions and water, and 
what is their effect? In particular we will 



focus on how the flow of ions and water can 
be controlled, i.e. we will investigate pos- 
sible gating mechanisms. This interest is 
based on the observation that the putative 
gate in many known ion channel structures 
(in particular, the nicotinic acetylcholine 
receptor nAChR [l3]jthe bacterial potas- 
sium channels KcsA jsf and KirBacl.l P], 
and the mechanosensitive channels MscL 
3 and MscS 0) is formed by a constric- 
tion made from hydrophobic residues. For 
example, the closed state nAChR struc- 
ture displays an ion pathway which is still 
wide enough (radius R ~ 0.31 nm) to ad- 
mit three water molecules (radius of a wa- 
ter molecule r^j = 0.14 nm) or one potas- 
sium ion (rK+ = 0.133 nm) with half of its 
ffist hydration shell intact. It is somewhat 
surprising that a pore does not have to be 
completely physically occluded to prevent 
the flow of ions. On the other hand ions 
can readily move through the KcsA |^ 
selectivity ffiter although its radius is less 
than 0.15 nm. When an ion enters the 
ffiter it has to shed its hydration shell at 
a high cost in free energy — the solvation 
free energy for a potassium ion is about 
-320 kJmol"^ O- The ffiter is lined by 
backbone oxygen atoms which coordinate 
the potassium ion and substitute for its 
hydration shell, thus reducing the desol- 
vation barrier of +320 kJmol"^ to about 
-M2kJmoH H. The putative gates dif- 
fer from the selectivity ffiter in that they 
are lined by hydrophobic residues. For 
nAChR [13, [TJ and MscL [3, [3 it was 
already hypothesised that these residues 
cannot substitute for water molecules so 
that the energetic cost of desolvation pre- 
vents the passage of the ion. 

It appears that the geometry of a pore 
(its radius, length and shape) and the 
chemical character of the pore wall have a 
great influence on the permeation of ions 
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and water. In addition, local flexibility (i.e. 
fluctuations in the protein structure as op- 
posed to concerted larger scale motions) of 
the pore lining might play a role, as seen 
in simulations of K"*" permeation through 
KcsA H. 

It seems difficult to comprehend the dy- 
namical nature of transport phenomena 
from static crystal structures alone. Com- 
putational methods can be used to comple- 
ment the experimentally observed picture. 
In particular, classical molecular dynamics 
(MD) can be used to investigate the be- 
haviour of water or ionic solutions in the 
environment presented by a protein. This 
realistic environment, however, makes it 
difficult to disentangle the contributions 
of various pore properties. In order to be 
able to reduce the number of parameters 
simplified pore models can be designed to 
cap ture the characteristics in question 
I25I ]. We focus on the influence of geome- 
try, pore wall character and flexibility. In 
previous studies we investigated the be- 
haviour of water in these "nanopores" . In 
particular, we found that below a criti- 
cal radius liquid water becomes unstable 
in the pore and the pore is predominantly 
filled with water vapour j25|. We hypoth- 
esised that a pore environment which can- 
not sustain liquid water would also present 
a high energetic barrier to an ion j2l|. Al- 
though it seems plausible that absence of 
water would imply the absence of ions it 
has not been demonstrated previously. 

In this work we explicitly consider ions 
in model pores and explore dimensions, 
pore wall character, and flexibility in more 
detail than in our previous work. We also 
present a simple thermodynamic model 
based on surface energies that explains the 
observed behaviour of water in hydropho- 
bic pores. 



2 Methods and Theory 

2.1 Model 

The pore models were constructed as de- 
scribed previously jl^. Briefly, they con- 
sist of concentric rings of methane-like 
pseudo-atoms of van der Waals radius 
0.195 nm (flgure Q). These are held in 
their equilibrium position by harmonic 
constraints with spring constant ko = 
1000 kJmol~^nm~^. A pore consists of 
two mouth regions (length 0.4 nm, i.e. 
one layer of pseudo atoms, and radius 
1.0 nm) at either end of the constriction 
site of length L (varied between 0.4 nm 
and 2.0 nm) and radius R (varied between 
0.15 nm and 1.0 nm). The pore was em- 
bedded in a membrane mimetic, a slab 
of pseudo atoms held on a cubic lattice 
with unit cell length 0.39 nm with har- 
monic springs of strength ko- Pores with 
a polar surface were created by placing 
partial charges of ±0.38e on atomic sites 
0.2 nm apart. The resulting dipole mo- 
ments pointed parallel to the pore axis 
with a magnitude of 3.6 D each, which is 
comparable to the dipole moment of the 
peptide bond (ca. 3.7 D jl^). 

2.2 Molecular Dynamics 

MD simulations were performed with 
GROMACS v3.1.4 and the SPC wa- 
ter model jl^]- The Lennard- Jones- 
parameters for the interaction between a 
methane-like pseudo-atom and the water 
oxygen are eco = 0.906493 kJ mol~^ 
and (Tco = 0.342692 nm; parameters for 
sodium and chloride ions are taken from 
the Gromacs force fleld. The integration 
time step was 2 fs and coordinates were 
saved every 2 ps. With periodic bound- 
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Figure 1: Simulation system: The solvent-accessible surface of a model pore with pore radius 
R = 0.55 nm is shown in light blue (parts of the pore close to the observer are removed for 
clarity). The pore is not completely hydrophobic due to two dipoles (red/blue) in the pore wall 
which have the same magnitude as the peptide bond dipole. Water molecules in the pore are 
depicted in licorice representation and as lines in the mouth region or the bulk (image created 
with VMD Q and Raster 3D Q). 



ary conditions, long range electrostatic in- 
teractions were computed with a parti- 
cle mesh Ewald method (real space cutoff 
1 nm, grid spacing 0.15 nm, 4th order in- 
terpolation j^) while the short range van 
der Waals forces were calculated within a 
radius of 1 nm. This treatment of electro- 
static effects is known to influence the per- 
meation of ions, especially in narrow pores. 
These artifacts need to be accounted for 
in order to obtain quantitative predictions 
of ionic currents S^]. However, though 
most likely present in our case, too, they 
are less important because of our emphasis 
on principles and order-of-magnitude ef- 
fects. The neighbour list (radius 1 nm) 



was updated every 10 steps. Weak cou- 
pling algorithms S^l were used to simulate 
at constant temperature (T = 300 K, time 
constant 0.1 ps) and pressure (P = 1 bar, 
compressibility 4.5 x 10~^ bar~^, time con- 
stant 1 ps) with the x and y dimensions 
of the simulation cell held fixed at 4 nm. 
The total thickness of the water reser- 
voir in the 2;-direction was 3.0 nm, en- 
suring bulk-like water behaviour far from 
the membrane mimetic. The initial sys- 
tem configuration was created by solvating 
the slab-embedded pore model with wa- 
ter. For simulations with ions, some wa- 
ter molecules were replaced with ions to 
reach the target concentration. Initially 
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there were always ions present in pores 
with R > 0.4 nm. A typical simulation 
box measured 3. 9x3. 9x4. 6 nm^ and con- 
tained about 1500 water molecules and 280 
pseudo atoms (and between 25 to 28 Na"*" 
and Cl~ ions each). 

2.3 State based analysis 

We define discrete states by mapping equi- 
librium states of the whole system (snap- 
shots from the MD equilibrium trajectory) 
onto numbers Wj. The same approach is 
used in statistical mechanics to map differ- 
ent microscopic states to one macroscopic 
state, characterised by the value of a state 
variable, here called u. This approach is 
used to label the two phase states that the 
pore water exhibits in our simulations. We 
either find liquid-filled pores or vapour- 
filled ones. Due to the small pore volumes 
"vapour" typically refers to zero or one wa- 
ter molecule in the cavity. We use the den- 
sity in the pore as an indicator of the phase 
state. As described previously [i^ we as- 
sign the state using a Schmitt-trigger pro- 
cedure ji^l in order to avoid spurious state 
changes due to fiuctuations from the inter- 
facial region. When the water density n{t) 
rises above 0.65 of the density of bulk wa- 
ter, no {no = 1.0 gcm~^ at T = 300 K and 
P = 1 bar) the liquid state, i.e. u = 1, 
is assigned to the phase state at time t. 
When n{t) drops below 0.25no the vapour 
state {uj = 0) is assigned. The pore in 
the liquid-filled state is termed "open" be- 
cause our simulations show that significant 
amounts of water pass through it; further- 
more, we also demonstrate in this work 
that a pore that sustains a liquid envi- 
ronment potentially allows ions to perme- 
ate. A vapour-plugged pore, however, will 
prevent ion permeation and is said to be 
"closed". 



Liquid-filled pores and vapour-filled 
pores are assumed to be in equilibrium. 
The "openness" (or probability for the oc- 
currence of the liquid, i.e. open, state is 

i^) = ^ dtuj{t) = -^, (1) 

sim Jo sim 

where Tsim denotes the total simulation 
time. To the total time that the pore is 
open, whereas Tc is the total time in the 
vapour or closed state. If (u) > ^ then 
equilibrium is on the side of liquid (and 
vapour is a meta stable state), otherwise 
the stable phase state is vapour. The equi- 
librium is governed by the equilibrium con- 
stant 

X(R) = '^"^^^ = ^sim - To{R) 

^ ' To{R) To{R) (2) 

= HR)r' - 1 

which is trivially related to the openness. 
[Note that for Tq = or Tc = it is mean- 
ingless to compute K because it indicates 
that the simulation time was too short to 
sample any state changes.] 

2.4 Thermodynamic model for 
liquid-vapour equilibrium in 
pores 

We investigate a subsystem of the whole 
simulation system comprising of the pore 
of volume V = LnR^. The subsystem can 
exchange water molecules with the bulk 
water outside the pore, which acts 
particle reservoir at average chemical po- 
tential /X. /X is implicitly determined by 
the constant average density of water in 
the bulk system, which in turn is pre- 
scribed by simulating at constant pres- 
sure. Hence we can use the grand poten- 
tial f2(T, V, yu) = —pV as the free energy to 
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describe the subsystem in the closed and 
open state 



24 



The equihbrium constant 
K is determined by the free energy differ- 
ence between the closed (vapour) and open 
(hquid) state, 



\xiK{R) 



(3) 



pAn{R) = (3 [n^{R) - QoiR)] = 
= -ln(Mi?))-^-l), 



using P = l/k^T. Following Evans [3f 
we write the free energy difference between 
the vapour and the liquid state with the 
corresponding surface contributions as 

An{R) = n^{R) - ni{R) 

= -p., (T, /i) L TiR^ + 2tiR L + 2t:R%, 
- {-pi{T,fi)L7iR^ + 27iRL-fi^) (4) 



Here the index w indicates the solid pore 
wall, I the liquid and v the vapour phase; 
for instance 7;^, is the surface tension or 
surface free energy per area of the liquid- 
vapour interface. The system is fairly close 
to bulk phase coexistence so we can ex- 
pand the pressure around the saturation 
chemical potential fisat{T) in a Taylor se- 
ries 



+ 



(5) 



A'=Msat 



where at saturation pi(T,fisa,t) = 
Pv(T, fisat)- This leads to a simple 
parabolic form for the free energy differ- 
ence between the two states: 
I 



An{R) = 2'ji^ - (/i - /isat) {nviT, /isat) - ni{T, /isat))^ 7ri?2 + 27rL(7„^ - 'ji^)R 

In. (6) 



7,^ + -AfiAn^iL 



TiR^ + 2itL A7^ R, 



where we define the distance of the state 
from saturation Afi := /i — /isat, the differ- 
ence in densities Auyi := ni — at sat- 
uration, and the difference in surface free 
energies of the two phases with the wall, 
^Iw ■= Ivw - liw The term AfiAn^i L 
is small for L < 10 nm as the system 
is close to phase coexistence (for L ^ 
1 nm it is about 10~^ times smaller than 
7i„ = 17 fc^Tnm"^ when estimated from 
An^iAuL ^ APL ^ 1 bar X L = 2.4 X 
lO-^fcfiTnm-^ xL,T = 300 K) and will be 
neglected. Only the difference between the 
surface free energies enters the model so we 
express it as the contact angle 9e, using the 
macroscopic definition from Young's equa- 
tion 'jvw—liw = liv COS 9e (see, for instance. 



Ref. j36|). Then equation IHl becomes 

Afi(i?,L,^e) = 2iTR-fi^{R + L cos^e)- 

(7) 

(equation 13 is similar to the simple model 
derived by Allen et al. [i^ but it in- 
cludes 'jvw — 'jiw and hence 9e instead of 
just 7;^.) For fixed pore length L and a 
given pore material, characterised by 6^, 
the graph of AQ{R) over the pore radius 
R describes a parabola containing the ori- 
gin. Although it is not a priori obvious 
that such a macroscopic treatment (equa- 
tion IHl or equation [7j) of a nanoscale system 
is meaningful the MD results for the free 
energy difference presented in section 13.11 
can be adequately explained within this 
model. 
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The connection between the model free 
energy AQ{R) and the behaviour of the 
system as observed in MD simulations, i.e. 
the openness {uj{R)), is established by in- 
verting equation ini yielding 

<^<«» = l+expMAfi(fi)] - 

For linear AQ{R) this represents a sig- 
moidal curve but non-linear terms (as in 
equation in} change its shape considerably. 

3 Results and Discussion 

First we investigate the behaviour of pure 
water in nanopores and quantify the influ- 
ence of geometry, pore surface and local 
flexibility on the equilibrium between liq- 
uid and vapour filled pores. Then simu- 
lations for a NaCl electrolyte are analysed 
and compared to the pure water case. The 
intention behind these studies is to explore 
more fully the ways in which a hydropho- 
bic gate may be opened, and how other 
factors might modulate such hydrophobic 
gating. This is particularly important in 
the context of current models for the gat- 
ing of nAChR [l^ that suggest that 
the transition between the closed and open 
states of the pore involves both an increase 
in pore radius and in polarity. We also 
note that recent discussion of the state 
(open vs. closed) of MscS have suggested a 
hydrophobic gating mechanism ISSll . simi- 
lar to that proposed for nAChR j2]J. 

3.1 Pure water 

A large number of simulations were run 
for water in model pores (total simula- 
tion time > 5 fis) in order to investigate 
the influence of pore radius (0.15 nm < 



R < 1 nm) and pore surface charac- 
ter (hydrophobic vs. amphipathic vs. hy- 
drophilic). Local flexibility (wall atom 
positional root mean square deviations 
(RMSDs) 0.039 nm < p < 0.192 nm) 
and temperature dependence in the range 
273 K < T < 450 K were examined for a 
hydrophobic R = 0.55 nm pore. 

Pore dimensions and surface character 

The water density in short nanopores os- 
cillates between liquid and vapour on a 
nanosecond time scale as seen in figure |21 
a manifestation of capillary evaporation 
and condensation at the nanoscale [25l |. 
In figure El the openness and the free en- 
ergy difference Afl between vapour and 
liquid states is shown for different pore 
surfaces. For hydrophobic and amphi- 
pathic pores a strong dependence of the 
pore state on the radius is apparent. The 
stable thermodynamic state switches from 
vapour {AQ < or {u) < i) to liquid 
[AQ > or (u;) > |) at a critical ra- 
dius Rc = —Lcos9e (using AQ{Rc) = 
in equation Ej). Rc is 0.56 ± 0.03 nm for 
the hydrophobic pore and 0.29 ± 0.02 nm 
for the amphipathic one. The functional 
form equation ini fits the data from the MD 
simulations (table well (the continuous 
lines in figure E)). The coefficient of the 
quadratic term, 7^ + ^Aft An^i L, is pos- 
itive and similar for both the hydropho- 
bic and the amphipathic (two dipoles, ab- 
breviated D2) pores (table [l]), consistent 
with the model equation IHl which pre- 
dicts this coefficient to be independent of 
the pore wall. For the polar pore (four 
dipoles, D4 ) the data is more ambiguous. 
Only if the first two data points at small 
radii are excluded from the fit as outliers 
is this coefficient positive (but still three 
times smaller than the one for the less hy- 
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60 80 100 
time t [ns] 

Figure 2: Liquid-vapour oscillations of water in a R = 0.65 nm apolar pore (top panel) when 
bathed in a 1.3 M NaCl solution. As indicated by the number N of ions within the pore, 
sodium (middle) and chloride ions (bottom) are only observed in the pore when there is also 
liquid water (n/nbuik ~ 0.8) present. Permeation events are indicated by triangles; ions do not 
permeate the pore during the vapour phases. 




Figure 3: Water in model pores. Left: openness (uj) over radius. Right: free energy difference 
Ail (in ksT) between vapour and liquid state. The grey region indicates radii smaller than 
the radius of a water molecule (0.14 nm). Data points are obtained from MD simulations 
with the errors estimated from block averages [3fll | . The continuous lines are fits of the model 
(equation EI) to the data points (right) or the openness computed from the model (equation |HI) . 
The vertical line indicates the radius of the closed nAChR gate {R = 0.31 nm). 
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Table 1: Parameters of the thermodynamic model equation 1^1 fitted to the MD results and 
resulting contact angle Be- Experimental values for 7/^, = 72 x 10"'^ Jm~^ = 17A;sTnm~^ 
and 9(, ~ 118° for water on a flat methyl (-CH3) terminated self assembled monolayer 
indicate that the simple thermodynamic model gives the right order of magnitude results. (All 
values at room temperature) 



pore surface character 












[k bT mn~'^] 


n 


bare CH4 "hydrophobic" 


+10±1 


-7.2±0.2 


134±3 


two dipoles D2 "amphipathic" 


+8.2±0.5 


-3.0±0.2 


111±2 


four dipoles "hydrophilic" 


+2.8±0.5 


+ 1.8±0.2 


51±10 



drophilic pores). The outliers show that 
a high density of local charges leads to a 
higher probability of the pore being liquid- 
filled than predicted by the macroscopic 
model, possibly indicating a shortcoming 
of the model to subsume inhomogeneous 
potentials into the surface tension terms. 
For our system parameters, the coefficient 
is in fact dominated by the water liquid- 
vapour surface tension 7^^ (see section 
and hence we will use its value as an ap- 
proximation to 7z„. A7^, the difference 
in surface tensions between the wall and 
vapour or liquid, becomes more positive 
with increasing polarity of the pore wall. 
It effectively measures the hydrophobicity 
of the wall. This becomes even more ap- 
parent when the (macroscopic) contact an- 
gle cos 6*6 = ^''^w/liv is formally computed 
(table P). Macroscopically, a hydropho- 
bic surface can be defined as one with 
9e > 90°. This allows us to call the ap- 
olar pore (6'e = 134°) "hydrophobic" com- 
pared to the "amphipathic" pore D2 {9^ = 
111°; still hydrophobic but with some "hy- 
drophilic" patches). A "hydrophilic" pore 
like the D4 system {6^ ~ 51°) is charac- 
terised by A'jm > or 6'e < 90° and liquid 
is always the preferred phase in the pore, 
regardless of R. 



Experimental macroscopic contact an- 
gles Be for water on fiat methyl (-CH3) ter- 
minated self assembled monolayers are re- 



411. A 



microscopic 



ported up to 118° |4C 
contact angle" of 135°±15° was calculated 
for a droplet of 90 water molecules (ra- 
dius R ~ 0.75 nm) on a methyl terminated 
film from MD simulations j42], consistent 
with our result of 6^ = 134°. MD sim- 
ulations seem to overestimate the experi- 
mental equilibrium contact angle. In this 
context we note that although the macro- 
scopic equilibrium contact angle is only es- 
tablished about 10 nm from the contact 
line, the local contact angle 61 in the core 
structure of the contact line (in the region 
3-10 nm from the contact line) in water 
(and electrolytes) is predicted to be larger 
than 9^ due to electrostatic screening ef- 
fects j36|. 

From the model free energy the open- 
ness can be calculated (equation IH} and 
displayed together with the MD results 
(figure El). Because the model effectively 
treats the liquid structureless contin- 
uum it is not meaningful to extrapolate to 
radii smaller than the radius of a water 
molecule (the shaded region in figure EI) . 
Nevertheless, the MD results lie on the 
model curve even in the direct vicinity of 
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this region. Further tests of equation [3 
by varying the length of the pore together 
with the radius confirm the model qualita- 
tively (data not shown). 

Our model implies that for both nano- 
and mesoscale pores the cost of creat- 
ing the liquid-vapour interface is the only 
force driving the filling of a hydropho- 
bic (A7^ < 0) pore. Little free energy 
A/iAn^i LtiR^ is gained by creating a bulk- 
like liquid in the pore instead of vapour. 



Local flexibility The influence of local 
fluctuations in protein structure can be 
modelled by changing the harmonic spring 
constant k that holds the pore wall atoms 
at their equilibrium positions. In a simpli- 
fied picture of a harmonic oscillator with 
an average energy ksT /2 per degree of 
freedom the RMSD p is directly related 



to k by p{k) = y (Armax)^ = \/^k^Tjk. 
This simple model overestimates the mea- 
sured RMSDs by only 10-20% (table 12), 
showing that the wall atoms behave like al- 
most independent, thermally driven, har- 
monic oscillators (figure E^.). At a fixed 
radius, increasing the flexibility (smaller 
/c, hence greater p) shifts the equilibrium 
towards vapour (figure . Conversely, a 
more rigid wall favours the condensation 
of water in the pore. Thus, water does 
not "push away" the pore walls to fill the 
pore but rather fluctuating pore walls ap- 
pear to disfavour the formation of adja- 
cent water layers. This effect cannot be 
explained with the assumption that the 
water molecules encounter a more nar- 
row pore "on average". Even if one as- 
sumes that the pore is narrowed by the 
RMSD of an atom down to an effective 
radius R{k) = Rq — p{k), the effect is 
still much stronger. When the thermo- 
dynamic model and the parameters for 



the hydrophobic pore (from tabled where 
/co = 4 X lO^/cfiT nm~^) are used to predict 
Af2(i?(/c)) then the prediction compared 
to the simulation results underestimates 
the effect for the more flexible pores but 
is in agreement for more rigid pores (fig- 
ure |3)d). We hypothesise that the random 
positional fluctuations of the wall atoms 
change the effective potential between the 
wall atom and water molecules by "smear- 
ing out" the attractive well of the inter- 
action potential. Preliminary calculations 
showed that an increase in flexibility (or 
temperature) would always decrease the 
well depth of the Lennard- Jones poten- 
tial, which is used to describe the inter- 
action between a water molecule and the 
methane molecules of the pore wall. Ni- 
jmeijer et al. 4^ demonstrated the depen- 
dence of the surface tension on the well 
depth e of the interaction potential analyt- 
ically. MD simulations of a water droplet 
with 1000 molecules {R ~ 2.5 nm) on 
a graphite-like surface with varying well 
depth ec-ow yielded a contact angle that 
varied from 30° to 180° 0. Thus it 
is the well depth of the interaction po- 
tential which ultimately determines the 
thermodynamic equilibrium in the pore. 
This is made explicit by considering equa- 
tion ini and noting that AQ{R) depends 
on the water- wall surface tensions through 
Ivw — liw, which in turn depend on e. For 
the case of water in confined geometries 
it was shown that capillary condensation 
and evaporation phenomena depend sensi- 
tively on e [2^ |4a] (or more precisely, 
on the effective well depth, the product of 
the wall atom density with e (25[). 

3.2 NaCI electrolyte 

As mentioned in the Introduction it has 
been hypothesised that a local hydropho- 
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I AQ.{R=0.55nm,k) [MD] 
AQ.{R,k) = Ail{R-Ap(k),kg) 




- 



G 
< 

CO. 



p(H=0.55nm,/c) [10' nm] 



Figure 4: Influence of flexibility on the liquid-vapour equilibrium, (a) For a hydrophobic 
R = 0.55 nm pore the harmonic restraint force constant k of the pore wall atoms was varied, 
resulting in a range of RMSDs p, which show the expected behaviour of thermic harmonic os- 
cillators. The free energy difference between liquid and vapour states of water in the pore 
shows strong dependence on the local flexibility p; rigid pores contain liquid water whereas 
flexible ones favour vapour. The effect can not be explained by a reduction of the effective 
pore radius by p alone [broken line in (b)]. All data at T = 300 K. 



bic environment would present a signifi- 
cant desolvation barrier to ion permeation. 
Our simulations of model pores bathed in 
a 1.3 M NaCl solution exhibit a striking 
change in the behaviour of the ions for 
pore radii R > 0.65 nm. Figure El shows 
that in pores of comparable dimensions to 
the closed nAChR pore {R = 0.35 nm) 
ions have a vanishing probability of enter- 
ing the pore, but at the open-state radius 
R = 0.65 nm the pore density rises to half 
the bulk value. The average density of ions 
in the centre of the pore (figure IH)) also ex- 
hibits a sharp increase near R = 0.65 nm, 
similar to the behaviour of the openness 
(figure El). 



The radially averaged densities show 
two pronounced water layers with an inter- 
layer distance d = 0.3 nm near the mem- 
brane mimetic. The ionic density drops 
off sharply at the maximum of the outer- 
most water layer and vanishes completely 
at the maximum of the water layer clos- 
est to the surface. Sodium and chloride 
ions are both present at the same distance 
from the surface. Because the membrane 
mimetic carries no charge and represents 
a low dielectric region (as do the less mo- 
bile bound water layers) it is energetically 
more favourable for the ions to be solvated 
in the high dielectric region of the bulk- 
like water. The bound water layers are 
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Figure 5: Equilibrium densities of water (top row), sodium ions (middle) and chloride ions 
(bottom) in hydrophobic model pores of radii R = 0.35 nm to 0.75 nm; the NaCl bulk con- 
centration is 1.3 moll~^. The narrow 0.35 nm pore mimics the closed gate of nAChR whereas 
the wide 0.65 nm pore approximates the open channel. Although wide enough to admit the 
ions physically, the narrow pore is effectively closed to ions ("hydrophobic gating"). The mem- 
brane mimetic is located in the horizontal region which appears black, i.e. void of water, in all 
images. 



12 




R [nm] 

Figure 6: Ionic density in the core region of hydrophobic pores. The core is taken to have a 
radius of about one ionic radius (0.1 nm) and exhibits the highest density of ions. The core 
density increases markedly at the critical radius R = 0.65 nm. The core density follows the 
openness (as shown for water in figure IH)). The dip in density at i? = 0.7 nm is due to very few 
state changes in the 150 ns simulation and a subsequent bias towards a 70 ns vapour period; 
in general the state switching frequency in the electrolyte systems is reduced compared to the 
pure water case (see also figure Bulk concentration of ions is 1.3 M. Lines are drawn as 
guides to the eye. 
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Table 2: RMSD of wall atoms in flexible pore models of radius R = 0.55 nm and influence on 
the liquid- vapour equilibrium of water in the pore at T = 300 K. 



restraint force constant k 


RMSD p{k) [10-^ nm] 


{uj) 


(3 An 


[10^ X fcsTnm-^] 


MD 


theory 


MD 


MD 


0.8 


1.7±0.3 


1.92 


O.OOliO.OOl 


-6.9±1.0 


2 


l.OiO.O 


1.22 


0.006±0.005 


-5.1±0.8 


4 


0.72±0.02 


0.86 


0.291±0.089 


-0.89±0.43 


8 


0.52±0.01 


0.61 


0.616±0.097 


0.47±0.41 


20 


0.35±0.01 


0.39 


0.755±0.060 


1.1±0.3 



not a consequence of the presence of the 
ions as they also form in pure water near 
a hydrophobic surface with the same in- 
terlayer spacing of 0.3 nm ji^, and in this 
case even a third water layer is discernible 
with d = 0.35 nm. Water ordering in 
cylindrical and slit pores in the absence 
of ions was also observed in simulations 
by other groups H H 113-0. Ef- 
fectively, ions are excluded from a zone of 
thickness 0.6 nm (Na"*") or 0.7 nm (Cl~) 
from the surface. The ion exclusion zone 
stretches into the pore without interrup- 
tion so that at R = 0.65 nm there is only 
a narrow channel of easily displaceable wa- 
ter which can admit ions. For wider pores, 
the ion density in the pore approaches bulk 
values (figure ini). 

Although the absence of ionic density 
can be taken ClS db strong indication for 
a closed pore the presence of equilibrium 
density does not necessarily imply a func- 
tionally open state, as the latter also re- 
quires rapid permeation of ions through 
a pore. As MD simulations output the 
trajectories of water and ions, the equi- 
librium flux ($o)) i-G- the total number 
of particles per nanosecond which com- 
pletely permeate the pore, can be mea- 
sured. Figure m shows $o for ions and wa- 
ter molecules (both for water in a 1.3 M 



NaCl electrolyte and for pure water) in 
hydrophobic pores of varying radius on 
a logarithmic scale. Water flux is about 
three orders of magnitude larger than ion 
flux and increases rapidly with increasing 
radii while R is smaller than the criti- 
cal radius Rc = 0.62 nm. Beyond Rc 
(when liquid becomes the stable phase) 
the slope in the logarithmic plot decreases 
from 15 nm~^ to 2.4 nm~^. The water cur- 
rent density ^q/{ttR'^) reaches a constant 
value of 300 ns~^ nm~^ so the increase for 
R > Rc is only due to the increase in 
pore diameter. The increase for R < Rc 
correlates with the openness (flgure EI), 
indicating that the main contribution to 
the flux stems from the open-state periods 
jisl . The presence of ions slightly shifts 
Rc from 0.56 ± 0.03 nm for pure water to 
0.62 ±0.07 nm for 1.3 M NaCl electrolyte. 
Qualitatively, this is explained by the in- 
crease of surface tension with ionic concen- 
tration 15(1. This leads to a higher contact 
angle 9e j^, which is also observed in our 
simulations where it increases from 134° to 
140°. Hence, Rc = —Lcos6e increases for 
electrolytes compared to pure water. 

The total number of ions that success- 
fully permeate the pore per nanosecond 
shows an increase over almost two orders 
of magnitude when the pore radius R is 
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R [nm] 

Figure 7: Equilibrium flux $0 of water and ions through hydrophobic model pores. The total 
number of ions that successfully permeate the pore per nanosecond shows an increase over 
almost two orders of magnitude when the pore radius R is increased from 0.6 nm to 0.65 nm. 
Water flux increases with the radius. The most notable influence of the 1.3 M NaCl electrolyte 
solution is the shift in the transition radius from 0.6 nm (only water) to 0.75 nm (water with 
NaCl). Lines are drawn to guide the eye and errors are estimated from block averages [s^. 



increased from 0.6 nm to 0.65 nm. Thus 
the pore "opens" at a radius much larger 
than the bare ionic radius. Ions perme- 
ate only when the pore is filled with liquid 
water (typical data in figure E)) and Na"*" 
permeation events are not correlated with 
Cl~ ions passing through the pore. 

Because the simulations are in thermo- 
dynamic equilibrium the net flux through 
the pore is zero for all R. In order to com- 
pare the MD equilibrium flux <l>0; which 
is determined by the intrinsic free energy 
barrier to ion permeation, to experimen- 
tally measured non-equilibrium fluxes $ of 
ions of charge g at a driving transmem- 
brane voltage V we employ rate theory and 
estimate 

«l>(V^;i?) = $o(i?)sinh^. (9) 
For R = 0.65 nm, V = 100 mV, q = le 



and T = 300 K we obtain $ ^ 1.2 ns~^ 
which is 40 times larger than the exper- 
imental value of 0.03 ns~^ for nAChR at 
0.2 M ionic concentration (calculated from 
a conductance of ca. 45 pS P). The dis- 
crepancy between our estimate and the ex- 
perimental data is not unexpected as equa- 
tion 121 has only qualitative character. Fur- 
thermore, nAChR presents a more com- 
plex pore lining surface than the model 
pore, and if the difference in bath concen- 
trations and access resistance was taken 
into account the discrepancy would likely 
be reduced by an order of magnitude. Nev- 
ertheless, the estimate demonstrates that 
the intrinsic barrier to ion permeation is 
small beyond the critical radius. 
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Figure 8: Temperature dependence of the liquid-vapour equilibrium in a hydrophobic 
R = 0.55 nm pore. Left: Change of the free energy landscape 0(n, T) with T. The density 
is calculated over the pore and the mouth regions which act as a buffer of constant chemical 
potential except for T = 273 K, where the mouth region also fills with vapour. Continuous 
lines are polynomial least square fits of order 4 (or order 6 for 300 K < T < 330 K), reminiscent 
of a Landau free energy with the water density in the pore and mouth region riMP as order 
parameter. The critical point of SPC water is near T = 373 K. Right: Openness and free 
energy difference between states Ail.. 



3.3 Sensing external parameters 

In the simplest model of sensing a sys- 
tem exists in an equilibrium with two 
states. External stimuli shift the equi- 
librium which elicits a cellular response. 
For example, a temperature-sensing chan- 
nel at normal temperature might be in the 
closed state. Elevated temperature shifts 
the equilibrium to the open state and the 
influx of ions initiates a signalling cascade 
that terminates in the sensation of heat. 
However, the exact mechanisms of gat- 
ing by temperature of e.g. TRP channels 
remain obscure js^]- We tested the in- 
fluence of temperature in the range from 
273 K to 450 K on water in the R = 
0.55 nm hydrophobic pore (data points 
for T > 373 K represent "supercritical" 



water because the critical point for the 
SPC water model is close to T = 373 K). 
The free energy landscape of the water- 
pore system depends on the temperature 
(figure IHk)- With increasing tempera- 
ture the liquid state (i.e. near water bulk 
density) decreases in stability. Accord- 
ingly, the openness decreases with temper- 
ature so that the vapour state, i.e. the 
closed stated, dominates at higher tem- 
peratures (figure IHb)- The decrease in 
openness with increasing T indicates that 
the difference of the wall surface tensions 
must become even more negative, 
i.e. vapour in contact with the wall is in- 
creasingly favoured. The effect would be 
similar to the one discussed for increased 
flexibility of wall atoms as the major ef- 
fect of temperature would be to smear 
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out the interaction potential between wa- 
ter molecules and wall atoms. This exam- 
ple demonstrates in principle how an ex- 
ternal signal, an increase in temperature 
from T = 300 K to 330 K, decreases the 
openness, and correspondingly the flux of 
water molecules, by a factor of six. 



4 Conclusions and outlook 

We have explicitly demonstrated a hy- 
drophobic gating mechanism for ions, us- 
ing simplified hydrophobic pores as mod- 
els for the closed gates of ion channels. 
There is a critical radius Rc above which 
a pore becomes effectively permeable to 
water or ions. Rc ~ 0.56 nm for water 
and Rc ~ 0.65 nm for ions, is much larger 
than the radius for a water molecule or 
a bare ion alone. This correlates nicely 
with the current view that the closed state 
structures of many ion channels contain 
gates formed by hydrophobic constriction 
sites. In nAChR the radius of the puta- 
tive gate is 0.31 nm, 0.17 nm in MscL, 
and 0.13 nm in KcsA. Models of KcsA and 
of Kv channels in their open states, based 
on the MthK and KvAP structures respec- 
tively, have gate radii of about 0.6 nm 5^, 
whereas the open-state nAChR structure 
opens up to about 0.65 nm. Our simula- 
tions show that the polarity of the pore 
wall can shift the critical radius consider- 
ably. This is reflected in the proposed gat- 
ing mechanism of nAChR ^3]. Not only 
does the radius of the constriction site in- 
crease but hydrophobic sidechains are also 
rotated out of the pore to expose the more 
polar, i.e. hydrophilic, peptide backbone. 
By combining hydrophobic gating with a 
change in surface polarity only a moder- 
ate change in radius is required to obtain 
a large physiological effect. Changes in lo- 



cal flexibility may also modulate the gat- 
ing behaviour. It is too early to ascertain 
their importance as such changes in flex- 
ibility have not been widely investigated 
experimentally for channels (or any other 
proteins). There are suggestions of possi- 
ble regulatory roles of changes in flexibility 
in e.g. KTN domains 5J], and in the mod- 
ulation of binding of TRAP protein to the 
Trp operon mRNA jH^. Temperature can 
also affect gating, in a fashion similar to 
local flexibility, but its direct effects seem 
too small to explain, for instance, the tem- 
perature sensitivity of TRP channels on 
their own. 

A thermodynamic model based on sur- 
face energies fits remarkably well the data 
of the atomic-scale MD simulation, even 
for very small radii. Though such a macro- 
scopic treatment is not a priori expected 
to give a satisfactory description of a mi- 
croscopic system there are other examples 
as for instance classical nucleation the- 
ory, which shares some similarities with 



our and others' models [2J,|56|. It gener- 
ally agrees quite well with experiments on 
the condensation of droplets from vapour, 
which implies that the use of macroscopic 
surface tension is valid even for droplets of 
radii of about 1 nm [53|. It appears that in 
this case water structure (i.e. a hydrogen 
bond network) is only important insofar 
that it is responsible for different wall-fluid 
surface tensions. Our results seem to cor- 
roborate the conclusions of Maibaum and 
Chandler jH^ that two-state behaviour of 
water in pores should only require a cold 
liquid close to phase coexistence and suf- 
ficiently different vapour-wall and liquid- 
wall surface tensions (which is, of course, 
where the "special" properties of water re- 
side in such a model). 

We have sketched out principles of gat- 
ing mechanisms in ion channels, based on 
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model channels. Ion channels in nature 
display complex conformational dynamics 
in relationship to gating. There is a need 
for better single molecule methods for ion 
channels to probe these phenomena exper- 
imentally. Combined electrical and op- 
tical methods look promising in this re- 
spect 5^ although considerable improve- 
ments in time resolution are still required. 
In the meantime there is a continued role 
for simulations and theory to enable us 
to bridge between static structure and dy- 
namic function. 

In addition to hydrophobic gating per 
se, our studies suggest the possible im- 
portance of gate flexibility in regulation 
of biological activity. The relationship 
between static structures, conformational 
change, and intrinsic flexibility in relation- 
ship to proteins and signalling merits fur- 
ther active investigation. In a recent re- 
view of protein-protein interactions and 
conformational changes, Goh et al. ET 



conclude that there is increasing support 
for a pre-existing equilibrium model. In 
such a model, proteins exist in a popula- 
tion of conformations, with ligand bind- 
ing leading to a change in the probabil- 
ity distribution of the ensemble. Some ev- 
idence in support of this have been ob- 
tained from recent simulation studies of 
ligand binding proteins (e.g. [6l[) but 
further (single-molecule) experimental and 
computational studies of a wider range of 
proteins are needed to more firmly estab- 
lish the general importance of changes in 
flexibility. 
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Abbreviations list 

MD molecular dynamics 

SPC simple point charge water model 

RMSD root mean square deviation 

nAChR nicotinic acetylcholine receptor 

MscS mechanosensitive channel of small con- 
ductance {Escherichia coli) 

MscL mechanosensitive channel of large con- 
ductance [Mycobacterium tuberculosis) 

KcsA potassium channel (Streptomyces livi- 
dans) 

KirBacl.l inward rectifier potassium chan- 
nel {Burkholderia pseudomallei) 

MthK calcium-gated potassium channel 
{Methanobacterium thermoautroph- 
icum) 

KvAP voltage-dependent potassium channel 
{Aeropyrum pernix) 



D2 model pore containing two dipoles (am- 
phipathic) 

D4 model pore containing four dipoles (hy- 
drophilic) 



TRP transient receptor potential channel 
family 



KTN potassium transport, nucleotide bind- 
ing domain 

TRAP trp RNA-binding attenuation protein 
trp tryptophan 
mRNA messenger RNA 
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Glossary 

Classical molecular dynamics A computa- 
tional scheme to calculate the time evo- 
lution of a microscopic system that is 
treated in atomic detail. Interactions 
are parametrised through classical two, 
three or four body potentials. A col- 
lection of these parameters is known as 
a "force field". The dynamics are ob- 
tained by integrating Newton's equa- 
tions of motion. 

Gating An ion channel can exist in an "open" 
state, which allows for ion permeation, 
and a "closed" state, when ions can- 
not pass the channel's gate and thus 
are prevented from crossing the mem- 
brane. Such a channel can be switched 
between states by external signals. 



Surface tension, surface free energy The 

free energy per area required to create 
an interface between two phases. It is 
always positive, i.e. it always costs free 
energy to create a surface compared to 
the homogeneous bulk phases. 



Contact angle The angle which is formed be- 
tween the tangent on a droplet of liq- 
uid where it is contact with the sur- 
face and the surface itself. The droplet 
is taken to be in equilibrium with its 
vapour. A surface material is called hy- 
drophobic ("water-hating") when the 
contact angle with water is greater than 
90° (the drop sits on the surface) and 
hydrophilic ( "water- loving" ) when it is 
smaller (the drop resembles a pancake) . 
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